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ABSTRACT 



The Naval Postgraduate School has an ongoing effort to study 
geomagnetic noise and micropulsations in the ULF frequency range 
(.05 4 f < 10 Hz). Data is collected by three orthogonally 
mounted coils at a remote land site and telemetered to the 
laboratory for computer analysis. 

To isolate data containing a micropulsation event, time 
series plots of the magnetic field were generated. The 
development of a double running average routine made possible 
the isolation of micropulsations in large data sets. A type 
PC-3 micropulsation was found and the coherence, ellipticity and 
polarization properties were determined as follows: Coherence = 

0.99, Degree of Polarization = 0.99 and the micropulsation was 
elliptically polarized. 

Power spectral density (PSD) plots summarizing data segments 
about two hours long occasionally contained structures found to 
be artificial and not representative of natural phenomena in the 
geomagnetic field. Methods to avoid this anomalous behavior in 
PSD plots are suggested. 
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INTRODUCTION 



I . 

This thesis is part of an ongoing effort at the Naval 
Postgraduate School to analyze ULF geomagnetic noise and 
micropulsations. The results of the Naval Postgraduate School 
studies could impact communications systems or systems in which 
geomagnetic noise and/or micropulsations introduce operational 
difficulties . 

An objective of this thesis was to develop the software 
necessary for obtaining computer generated plots of geomagnetic 
field versus time. Then, with the plots, sections of data 
containing obvious geomagnetic micropulsations were located. 
Those sections were analyzed to determine the power spectral 
density (PSD), coherence, degree of polarization and 
ellipticity of the micropulsation using already developed 
software [Ref. 1]. 

After analyzing some of the plots (PSD, coherence, etc.) 
produced by the software documented in [Ref. 1], it was 
demonstrated that problems existed in the software. New 
software was developed to obtain times series raw voltage and 
geomagnetic field plots. Validation of the software required 
an on-site experiment. Documented within this thesis is the 
software to produce time series voltage and geomagnetic field 
plots and the questions the validation experiment raised 
concerning the sensing system and software. 
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Using the newly developed software, a micropulsation event 
was located. A double running average software routine 
accentuated the micropulsation to permit the determination of 
its frequency and amplitude. The micropulsation event was 
further analyzed to obtain information on coherence, degree of 
polarization and ellipticity. 
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II . BACKGROUND 



A. COLLECTION SYSTEM 

The Naval Postgraduate School has two operating geomagnetic 
sensing systems. The land site is located at the La Mesa 
Village Housing Area and has three orthogonally mounted 
induction coils. The second operating system contains two 
orthogonally mounted induction coils that can be placed on the 
bottom of Monterey Bay. When both systems are operating, the 
data is collected simultaneously. The ground work has been laid 
for a third system located at Chew's Ridge. Three orthogonal 
and simultaneous measurements can be taken only between two 
sites since the School has only five induction coils and a sixth 
just calibrated. This thesis used data from the La Mesa Village 
site and the software associated with it. 

An induction coil senses the geomagnetic field fluctuations 
based on Faraday's Law of induction, induced emf =- » which 

describes the relation between the induced emf and the time 
rate of change of magnetic flux through the coil. N is the 
number of turns in the coil. 

The voltage induced in the coils is sampled 32 or 64 
times/second and is amplified approximately one million times. 
The voltage data undergoes pulse-code-modulation (PCM) for 
noiseless VHF radio link transmission from the La Mesa 
Village site, to the recording system in the Geophysics Signal 
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Processing Laboratory located in Spanagel 531. Data reduction 
on the Naval Postgraduate School's IBM 3033 computer follows 
digitization of the PCM data. A block diagram of the system is 
provided in Figure 2.1. A more detailed description of the 
system is given in Reference 2. 

1 • System Calibration 

For purposes of documentation, the calibration of the 
third coil for either the Chew's Ridge or Monterey Bay system 
is being included in this thesis. The general method of 
induction coil calibration will be covered here. For a more 
indepth explanation, see Reference 1. 

The sensing coil and its associated electronics are 
calibrated using a Helmholtz coil apparatus to establish a 
uniform magnetic field. By placing the sensing coil into the 
field produced by the Helmnolcz coil, one can establish a 
relationship between the Helmholtz coil magnetic field magnitude 
and frequency and the coil system voltage. Once the response of 
the coil system is known for the frequency band of interest, it 
can be included in the system transfer function used in data 
reduction . 

The experiment to determine the sixth coil system 
transfer function was done at the La Mesa Village site. A 
Wavetek Model 142 signal generator supplied a sinusoidal current 
to a 1.22 meter diameter, 0.61 meter high Helmholtz coil. The 
current was measured across a 996 ohm resistor in series with 
the Helmholtz coil. A Hewlett Packard HP-3582A spectrum 
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Figure 2.1. System Data Flow 



analyzer measured the output voltage of the coil system as a 
function of frequency. 

The experiment measured the response of the coil system 
in the frequency band .05 hertz to 20 hertz for applied fields 
of 0.02, 0.2 and 1.0 nanoteslas. The coil system response is 
shown in Figure 2.2. The transfer function of this coil system 
corresponds closely to those of coil systems previously 
calibrated. The only noticable difference is that the 0.02 
nanotesla response at 15 hertz is approximately 15 volts/ 
nanotesla higher than those measured before [Ref. 1]. It 
should be noted that this difference will not influence data 
analysis in the frequency range of interest. The transfer 
function algorithm is lis'ted in Table 2.1. 

B . PREVIOUS SOFTWARE 

Software has been developed (see References 1 and 3), to 
extract power spectral densities for each coil, coherence 
between the coils, degree of polarization in the three measure- 
ment planes, ellipticity in the three measurement planes and a 
variety of Stokes parameters. The following is a basic 
description of the computer code that produces the above 
mentioned plots. First, the voltage data is read off the 
digital tape using a subroutine called RD provided by Dr. Tim 
Stanton of the Oceanography Department. A parameter ISEC, 
representing the number of seconds one wishes to advance the 
tape, is frequently used. ISEC establishes the number of 
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Figure 2.2. Coil Transfer Function 



Table 2.1. System Transfer Function Algorithm 



DC 9 L=1 ,N 
rFQ=FREQ(L) 

IF(FRQ. LE. 25 . )G0 TO 1 
XX ( l )=XX (L ) / 28. 

GC TG 3 

1 Ir( FRQ . LE. 15 . )Gu TO 2 
XX(L)=XX (l ) / ( 105 .5-3.14*FRQ) 

VY(L) = YV (L) / (131 .32-7.536*FRQ) 

ZZ ( L )= ZZ (L ) / ( 177.26- 7. 484*FRQ) 

G C TO 8 

2 Ir( FRQ . LE. 10. )GO TO 3 
XX(L)=XX(l)/ ( 5. 9 53* FRQ- 30. 97) 

YY ( L) = YY (L) / (7. 166* FRQ -3 9.99 ) 

ZZ( L )= ZZ (L ) / ( 6.4 9* FRQ- 32 .35) 

C-C TO 3 

3 IF (FRQ. LE. 7. 5)GG TO 4 
XX(L)=XX(L) / (3.492* FftQ-6. 31 ) 

YY(L)=YY ( L ) / (4. 2 52* FRQ- 10. o5) 

ZZ( L|=ZZ (L; / (4.04** FRQ -7 .o9) 

V3 L « U 0 

I P I r" ft** • L O • /sjiJ 1 G 0 

XX ( L ) = XX ( L ) / (2.6311 *FR 3+0. 14oo7) 

YY(L)=YY(L)/(3.0 12* FRQ- 1 . 55) 

ZZ( L) = ZZ (L) / (3.134* FRQ-i .44 ) 

GC TC 3 

5 I F ( FRQ . L E. 3 . )GCi TO 6 

XX ( l) = XX(l ) / (2.6311 *FRQ+0. Woo 7) 

7 Y Y i L )= Y Y (L )/ (2.702*FRQ) 

ZZ( L) = zz (L) / (2.92*FRQ) 

GC TO 3 

o XX ( L )= XX (L)/ (2.72*FRQ) 

GC TC 7 
3 CONTINUE 

TF(L)=(XX(L )*C0S0 + YY (L ) *CuS01 ) *Cu SoO + ZZ(L)*CCS3G 
9 CCNTINJE 



15 



times RD is called and the data not stored. Previously, tape 

advances of approximately 30 seconds or ISEC = 30 were -used. 

Once the tape is advanced the desired amount, data analysis can 

begin. Each digital tape contains about 90 minutes or 1.728 x 
s 

10 pieces of data per axis resulting from the sampling rate of 
32 samples/second. Because of the large amount of data, the 
analysis has to be accomplished in blocks. Data analysis takes 
place on 256 seconds or 8192 frames of data at a time. RD is 
called 8192 times and places the voltages for each coil in an 
array. The voltages are integer values from 0 to 4096 
representing -5 volts to +5 volts . Ideally the value 2048 
represents 0 volts. The integers between 0-4096 are normalized 
to + 5 volts and placed in new complex arrays. These arrays 
containing 8192 pieces of data are Fourier transformed to the 
frequency domain where the system transfer function is applied. 
After applying the transfer function, the calculations begin 
for a determination of PSD's, coherence, ellipticity, etc. for 
the first block. After they are complete, the next block of 
8192 pieces of data starts through the program. Once the 
desired amount of data has been processed, a NONIMSL subroutine 
called DRAVJP plots the data. For a more complete description 
of this software, see Reference 1, Reference 3 and 
Reference 4. 
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C. ANALYSIS OF INITIAL DATA PLOTS 



Many of the PSD's, coherences, ellipticities and degree of 
polarization plots produced by the software of Reference 1 and 
Reference 3 have unexpected characteristics. The PSD's show 
“humps" that are shown in Figures 2.3-2. 5. These “humps" are 
not characteristic in "normal" PSD's. The coherence plots show 
a coherence of essentially one degrading to hash in the 
frequencies greater than 2 hertz. Figures 2. 6-2. 8. We see a 
similar behavior in the plots of degree of polarization. These 
plots indicate a very highly polarized field from 0.02 hertz to 
1.0 hertz. Figures 2.9-2.11. The ellipticity plots show a very 
linearly polarized field for these frequencies in Figures 
2.12-2.14. 

After examining the time series magnetic field plots, we 
noted a section of very large magnitude fluctuations (about one 
to three orders of magnitude greater than the general back- 
ground) at the beginning of the digital tapes Figures 2.15-2.17. 
We found that if the digital tape was advanced past these large 
fluctuations the output (PSD, coherence, etc.) appeared 
"normal". Figures 2.18-2.29. If the program was executed on 
just the section containing the large fluctuations, the PSD's 
had the "humped" characteristic, the coherences were one, the 
degree of polarization was very high and the ellipticity was 
linear. Figures 2.30-2.41. This demonstrates that these large 
fluctations dominate the output of an analysis in the frequency 
domain. A tape advance of 300 seconds or ISEC = 300 is 
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Figure 2.3. PSD X-Coil, 17 August 82, 2240-2348 Local. 

Amplitude in dB (REF nT**2/Hz) (20 units/in) vs. Log 

Frequency (Hz) (0.5 units/in), 16 Averages 
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Figure 2.5. PSD Z-Coil, 17 August 82, 2240-2348 Local. 

Amplitude in dB (REF nT**2/Hz) (20 units/in) vs 

Log Frequency (Hz) (0.5 units/in), 16 Averages 
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Figure 2-6* Coherence X and Y coils, 17 August 82, 2240-2348 Local. 

Coherence X and Y Coil (0.2 units/in) vs. Log Frequency 
(Hz) (0.5 units/in), 16 Averages 
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Figure 2.7. Coherence Y and Z Coils, 17 August 82, 2240-2348 Local. 

Coherence Y and Z Coils (0.2 units/in) vs. Log Frequency 
(Hz) (0.5 units/in), 16 Averages 
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Figure 2.8. Coherence Z and Z Coils, 17 August 82, 2240-2348 Local. 

Coherence Z and X Coils (0.2 units/in) vs. Log Frequency 
(Hz) (0.5 units/in), 16 Averages 
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Figure 2.9. Degree of Polarization X-Y Plane, 17 August 82, 2240-2348 Local. 

Degree of Polarization (0.2 units/in) vs. Log Frequency (Hz) 

(0.5 units/in), 16 Averages 




25 



Figure 2.10. Degree of Polarization Y-Z Plane, 17 August 82, 2240-2348 Local. 

Degree of Polar iz^ltion (0.2 units/in) vs. Log 
Frequency (Hz) (0.5 units/in), 16 Averages 
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Figure 2.11. Degree of Polarization Z-X Plane, 17 August 82, 2240-2348 Local. 
Degree of Polarization (0.2 units/in) vs. Log Frequency 
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Figure 2.13. Ellipticity Y-Z Plane, 17 August 82, 2240-2348 Local. 

Ellipticity (0.5 units/in) vs. Log Frequency (Hz) (0.5 units/in) 




29 



Figure 2.14. Ellipticity Z-X Plane, 17 August 82, 2240-2348 Local. 

Ellipticity (0.5 units/in) vs. Log Frequency (Hz) (0.5 units/in) 
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Figure 2.15. X-Coil Magnetic Field, 17 August 82, 2240-2348- Local. 

Amplitude (nanoteslas : 10 units/in) vs. Time (seconds : 500 units/in) 
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Figure 2.16. Y-Coil Magnetic Field, 17 August 82, 2240-2348 Local. 

Amplitude (nanoteslas : 10 units/in) vs. Time (seconds : 500 units/in) 
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Figure 2.17. Z-Coil Magnetic Field, 17 August 82, 2240-2348 Local. 

Amplitude (nanoteslas : 5 units/in) vs. Time (seconds : 500 units/in) 
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Figure 2.21. Coherence X-Y coils, 17 August 82, 2245-2353 Local. 

Coherence X-Y Coils (0.2 units/in) vs. Log Frequency (Hz) 
(0.5 units/in), 16 Averages 
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Figure 2.22. Coherence Y-Z Coils, 17 August 82, 2245-2353 Local 

Coherence X-Y Coils (0.2 units/in) vs. Log Frequency (Ilz) 
(0.5 units/in), 16 Averages 
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Figure 2.23. Coherence Z-X Coils, 17 August 82, 2245-2353 Local. 

Coherence Z-X Coils (0.2 units/in) vs. Log Frequency (Hz) 
(0.5 units/in), 16 Averages 





39 



Figure 2.24. Degree of Polarization X-Y Plane, 17 August 82, 2245-2353 Local. 

Degree of Polarization (0.2 units/in) vs. Log Frequency (Hz) 

(0.5 units/in), 16 Averages 
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Figure 2.25. Degree of Polarization Y-Z Plane, 17 August 82, 2245-2353 Local. 

Degree of Polarization (0.2 units/in) vs. Log Frequency (Hz) 

(0.5 units/in), 16 Averages 
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Figure 2.26. Degree of Polarization Z-X Plane, 17 August 82, 2245-2353 Local. 

Degree of Polarization (0.2 units/in) vs. Log Frequency (Hz) 

(0.5 units/in), 16 Averages 
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Figure 2.28. Ellipticity Y-Z Plane, 17 August 82, 2245-2353 Local. 

Elliptici ty (0.5 units/in) vs. Log Frequency (Hz) (0.5 units/in) 
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Figure 2.29. Ellipticity Z-X Plane, 17 August 82, 2245-2353 Local. 

Ellipticity (0.5 units/in) vs. Log Frequency (Hz) 

(0.5 units/in), 16 Averages 
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Figure 2.33. Coherence X-Y Coils, 17 August 82, 2240-2250 Local. 

Coherence X-Y Coils (0.2 units/in) vs. Log Frequency (Hz) 
(0.5 units/in), 2 Averages 
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Figure 2.34. Coherence Y-Z Coils, 17 August 82, 2240-2250 Local. 

Coherence Y-Z Coils (0.2 units/in) vs. Log Frequency (Hz) 
(0.5 units/in), 2 Averages 
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Figure 2.35. Coherence Z-X Coils, 17 August 82, 2240-2250 Local. 

Coherence Z-X Coils (0.2 units/in) vs. Log Frequency (Hz) 
(0.5 units/in), 2 Averages 
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Figure 2.36. Degree of Polarization X-Y Plane, 17 August 82, 2240-2250 Local. 

Degree of Polarization (0.2 units/in) vs. Log Frequency (IIz) 

(0.5 units/in), 2 Averages 
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Figure 2.37. Degree of Polarization Y-X Plane, 17 August 82, 2240-2250 Local 
Degree of Polarization (0.2 units/in) vs. Log Frequency (Hz) 
(0.5 units/in), 2 Averages 
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Figure 2.38. Degree of Polarization Z-X Plane, 17 August 82, 2240-2250 Local. 

Degree of Polarization (0.2 units/in) vs. Log Frequency (Hz) 

(0.5 units/in), 2 Averages 
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Figure 2.39. Ellipticity X-Y Plane, 17 August 82, 2240-2250 Local. 

Ellipticity (0.5 units/in) vs. Log Frequency (Hz) (0.5 units/in) 
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Ellipticity (0.5 units/in) vs. Log Frequency (Hz) (0.5 units/Ln) 




56 



Elliptici ty (0.5 units/in) vs. Log Frequency (Hz) 



recommended at the beginning of the programs producing the above 
mentioned plots. Time series plots should also be made to 
provide a further check. 

An explanation of the source of the large fluctuations is 
that the PCM to digital conversion begins before the actual PCM 
data starts on the tape. The decoding equipment considers the 
noise on the section of tape before the PCM data starts as 
signal . 

The above analysis clearly proves that the "hump" like 
structure in the PSD's and the "nice” behavior observed in the 
corresponding coherence, ellipticity and degree of polarization 
are artifical and do not represent any real phenomena in the 
geomagnetic field. Thus, the data in Figures 2.3-2.14 and 
Figures 2.30-2.41 does not appear to be reliable. Unfortu- 
nately, this artificially induced, and therefore unphysical, 
behavior of the polarization has been published recently 
[Ref. 7]. The conclusion found in Reference 7 that 

Geomagnetic fluctuations in land are well polarized 
below 1.0 Hz and that the polarization is quite 
frequently, though not consistently, linear. The 
situation under the sea is similar. Furthermore, 
a land-sea coherence study in the horizontal plane 
showed very high coherence below 1.0 Hz. 

are not substantiated by the physical data. Plots of the 

artificially induced behavior are published in References 1, 2, 

and 3 . 
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The 4 Hz peak present in the PSD's shown in Figures 2.13- 
2.20 is the result of 60 Hz aliasing. 60 Hz magnetic fields 
generated by local power lines exist at the measurement site. 
Digitization at a 32 Hz or 64 Hz rate of the analog voltages 
from the sensing coils shifts the true frequency of the power 
line signals down at 4 Hz. The 60 Hz aliasing results in the 
anomalis behavior at 4 Hz seen in the coherence, degree of 
polarization and ellipticity plots shown in Figures 2.21-2.29. 

Another peak is present at about 1.2 Hz for the PSD's in 
Figures 2.18-2.20. For this data set, the 1.2 Hz signal 
correlates with a 1.2 Hz signal identified in the time series 
data. The origin of the 1.2 Hz signal is unknown. The 4 Hz and 
1.2 Hz peaks are of little concern in this study since both are 
well above the frequency of the micropulsation analyzed. 

Another "suspicious" feature of the PSD's is the large 
amount of hash appearing from 2 to 10 hertz. We will gain 
further insight into this problem in analyzing the time series 
voltage data. 

From the analysis of the abnormal behavior in the PSD's, 
coherence, etc., future coil data sets should receive a 
balanced treatment between the time domain and the frequency 
domain to avoid misinterpretation. 
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III. TIME SERIES DATA 



A. TIME SERIES VOLTAGE DATA 
1 . Voltage Software 

The importance of time series data as a check became 
evident when we uncovered the problem with the initial PSD's, 
coherences, etc. Computer code was developed to produce both 
unsmoothed and smoothed times series voltage plots. The reason 
for producing both smoothed and unsmoothed voltage plots will 
become apparent in the time series voltage analysis section. 
Time series voltage plots could provide a means of monitoring 
the data for system "glitches” which might later produce 
erroneous ?SD‘s, etc., and suspicous features on time series 
magnetic field plots. 

The Fortran computer program that produces the voltage 
plots is compiled and run on the Naval Postgraduate School IBM 
3033 VM or "batch" system. Both the unsmoothed voltage program 
(VOLTR ) and the smoothed voltage program (VOLTS) require 2.048 
megabytes of core memory and approximately 15 minutes of 
central processor unit (CPU) time to run. These figures 
provide a comfortable margin to insure that the program does 
not "bomb" due to lack of core or run time. 

The design of both VOLTR and VOLTS is the same with the 
exception that a smoothing algorithm is applied to the data in 
the case of VOLTS. Both programs have incorporated in them the 
digital tape advance package using ISEC as the number of 
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seconds one wishes to advance the tape. Data is read off the 
digital tape in blocks of 8192 frames by the subroutine RD 
referenced earlier. As before, this data is in integer form 
and represents voltages between 4^5 volts by numbers between 0 
and 4096. The integer representing 0 volts is ideally 2048, 
however, I found that the "real" zero values wander. The 
reason for this "wander" is not understood. For good results, 
the following values worked well for the respective coils. 



X-coil : 


zero = 


1966 


Y-coil : 


zero = 


2085 


Z-coil : 


zero = 


2539 



The voltage values normalized to +5 volts are placed in new 
real arrays. At this point the next block of data is read off 
the digital tape and goes through the same process. VOLTR and 
VOLTS analyze eight blocks or 34 minutes of data. The linear 
arrays holding the voltage data are dimensioned to 65,536 
elements given a sampling rate of 32 samples/second. For every 
one of the 65,536 voltage values there is a time value in 
another array. Thus four arrays are required to handle the 
data, three arrays containing voltage data and one that 
contains the respective times for each of the measurements. At 
this point, the VOLTR program calls the NONIMSL subroutine 
DRAWP to plot the data. The VOLTS program uses a smoothing 
algorithm to reduce the amount of high frequency hash. 

The smoothing algorithm is a double running point 
average executed on all three axes of data. Since I was 
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interested in periods of approximately 10-45 seconds, I 
averaged over 144 points. With a running average over 23 
percent of the oscillation, the unwanted background can be 
removed without destroying the oscillations. Averages greater 
than 50 percent of the period would partially remove the 
oscillation of interest. The construction of the algorithm 
makes changing the number of points averaged over trivial. 

One need only change a loop index and three divisors to obtain 
a different average. The smoothing algorithm is provided as 
Table 3.1. The computer code for both VOLTR and VOLTS is 
attached in Appendix A and Appendix B. 

2. Time Series Voltage Data Analysis 

One can gain a great deal of insight into the contents 
of the data by looking at the smoothed and unsmoothed voltage 
products. Figures 3. 1-3. 3 are unsmoothed voltage plots 
representing the first 34 minutes of the digital data tape GMDT 
11 recorded 17 August 82, 1301-1408 local time. A one minute 

tape advance was used in producing this data. Immediately one 
notices the large amount of 60 hertz aliasing in the time series 
unsmoothed voltage. The 60 hertz aliasing appears as 4 hertz 
because of the 32 sample/second sampling rate. The sensing 
system has a low pass filter that should be improved to remove 
the aliasing problem. The 4 hertz background can be seen more 
clearly in Figures 3.4— 3.6. Above the 60 hertz aliasing 
background, we see voltage spikes ranging in magnitude from 1 
volt to over 4 volts. Expanding the time scale to 50 
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Table 3.1. Smoothing Algorithm 



CC 73 L2 =1 » 2 

- -Q 

DC 7<t IS=1, 65318 
SUMX=0 . 0 
SUM Y=Q . 0 
SUM2=0 . 0 
CG 75 J= 1» 144 
SUMX=ZZXI(Q+J1+SUMX 
SLMY=ZZY1(Q+ JI + SUMY 
SUMZ=ZZ V 1( Q+ J I+SUMZ 
75 CCNTINUE 

ZZXKI S) =SUMX/144. 
ZZY 1(1 S )=SUMY/l44. 
ZZVKI SJ -=$UM l/l+H, 
*=G + 1 

74 CCNTINUE 
73 CCNTINUE 
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Figure 3.1. X Coil Voltage, 17 August 82, 1302-1336 Local. 

Amplitude (volts : 2 units/in) vs. Time (seconds : 500 units/in) 
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Figure 3.2. Y Coil Voltage, 17 August 82, 1302-1336 Local. 

Amplitude (volts : 2 units/in) vs. Time (seconds : 500 units/in) 
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Figure 3.3. Z Coil Voltage, 17 August 82, 1302-1336 Local. 

Amplitude (volts : 2 units/in) vs. Time (seconds : 500 units/in) 
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Figure 3.4. X-Coil Voltage, 17 August 82, 1322-1323 Local. 

Amplitude (Volts: 2 units/in) vs. Time (seconds : 10 units/in) 
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Figure 3.5. Y-Coil Voltage, 17 August 82, 1322-1323 Local. 

Amplitude (volts : 2 units/in) vs. Time (seconds : 10 units/in) 
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Figure 3.6. Z-Coil Voltage, 17 August 82, 1322-1323 Local. 

Amplitude (volts : 2 units/in) vs. Time (seconds : 10 units/in) 



seconds of data, the structure of the voltage spikes becomes 
apparent in Figures 3. 4-3. 6. The "square well" shape is 
characteristic of PCM dropouts or instances during the PCM to 
digital conversion synchronization is lost between the tape 
recorder containing the PCM tape and the decoding unit. We do 
not know whether PCM dropouts always occur on all three 
channels, but current evidence suggests dropouts are 
simultanoues on all channels. 

Figures 3.7-3. 9 again show the first 34 minutes of GMDT 
11 only smoothed by the 144 point running average. The smooth- 
ing algorithm removes the 60 hertz alising. The PCM dropouts 
appear as "glitches" in the otherwise quasi sinusoid charac- 
teristic of a micropulsation. 

We should note that the unusually large amount of hash 
present between 2 hertz and 10 hertz in the PSD 1 s and the 
average slope of the PSD curve for those frequencies may in part 
be due to the PCM dropouts that occur on the tapes . A 
"deglitching" algorithm could be developed to remove the sudden 
large amplitude structures injected by the system into the 
voltage data if hardware modifications fail. 

B. TIME SERIES MAGNETIC FIELD DATA 
1. Magnetic Field Software 

The time series magnetic field software is basically 
the same as that which produces the PSD's, etc. [Ref. 1 and 
Ref. 3]. The differences begin after the forward Fourier 
transform and the application of the transfer function. We 
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Figure 3.7. X Coil Voltage, 17 August 82, 1302-1336 Local. 

Amplitude (volts : 0.2 units/in) vs. Time (seconds : 500 units/in 
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Figure 3.8. Y Coil Voltage, 17 August 82, 1302-1336 Local. 

Anplitude (volts : 0.1 units/in) vs. Time (seconds : 500 units/in) 
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Figure 3.9. Z Coil Voltage, 17 August 82, 1302-1336 Local. 

Amplitude (volts : 0.1 units/in) vs. Time (seconds : 500 units/in) 
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Figure 3.10a. X-Coil Magnetic Field, 17 August 82, 1302-1336 Local. 

Amplitude (nanoteslas : 5 units/in) vs. Time (seconds : 500 units/in) 
8 point smoothing 
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Figure 3.12a. Z-Coil Magnetic Field, 17 August 82, 1302-1336 Local. 

Amplitude (nanoteslas : 2 units/in) vs. Time (seconds : 500 uriits/in) 
8 point smoothing 
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Figure 3.10b. X Coil Magnetic Field, 17 August 82, 1302-1336 Local. 

Amplitude (nanoteslas : 5 units/in) vs. Time (seconds : 500 units/in) 
144 point smoothing 
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Figure 3.11b. Y Coil Magnetic Field, 17 August 82, 1302-1336 Local. 

Amplitude (nanoteslas : 1 unit/in) vs. Time (seconds : 500 units/in) 
144 point smoothing 
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Figure 3.12b. Z Coil Magnetic Field, 17 August 82, 1302-1336 Local, 

Amplitude (nanoteslas : 1 unit/in) vs. Time (seconds : 500 units/in) 
144 point smoothing 



should note that the transfer function as designed ignores 
phase information since it is just a least squares fit- to data 
like that of Figure 2.2. Once the transfer function is 
applied, the reverse transform is performed. We obtain the 
complex absolute values of the magnetic field data (in 
nanoteslas) and place the data in new real arrays. To insure 
that there are no discontinuities between blocks of data, 

Figures 3.10a, 3.11a and 3.12a, all the blocks are connected to 

the first by adding or subtracting a constant to each element, 
Figures 3.10b, 3.11b, 3.12b. The data now undergoes the same 

144 point double running average as the voltage to bring out 
fluctuations with periods of 10-45 seconds. Again DRAWP is used 
to plot the data. The computer program LVFTC1 produces the 
smoothed magnetic field plots. It is provided in 
Appendix C. 

2 . Smoothed Magnetic Field Data Analysis 

Again looking at the first 34 minutes of digital data 
tape GMDT 11, we see in the magnetic field data the same 
exceptional structures as seen in the voltage data, Figures 
3.10-3.12. The glitches, however, do highlight one problem 
rather well. From Figures 3. 7-3.9 we note that positive 
voltage spikes produce negative magnetic field spikes. This 
would indicate a pi phase shift is introduced either in the 
transfer functions or in the digital Fourier transform. After 
investigating another digital tape GMDT 3A, 18 August 82, 

0121-0251 local, we see another behavior. Figures 3.13-3.15 
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Figure 3.13. X Coil Voltage, 17 August 82, 2241-2315 Local. 

Amplitude (volts : 0.1 units/in) vs. Time (seconds : 500 units/in) 
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Figure 3.14. Y Coil Voltage, 17 Voltage 82, 2241-2315 Local. 

Amplitude (volts : 0.1 units/in) vs. Time (seconds : 500 units/in) 
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are the smoothed voltages for the first 34 minutes of data 
on GMDT3A after a one minute tape advance. Figures 3.-16-3.18 
are the smoothed magnetic fields for the same time period. The 
X-coil voltages show three positive spikes while the X-coil 
magnetic field shows three negative spikes. This supports the 
theory of a pi phase shift. However, the Y-coil voltage shows 
two positive spikes and one negative spike, while the Y-coil 
magnetic field has three positive spikes. Thus it appears we 
have a rather arbitrary phase problem that may find its origin 
in the transfer functions that do not maintain the correct 
phase . 

From Figures 3.16-3.18 we note exceptional features 
that do not correspond to voltage spikes. The exceptional 
features are separated by time periods of 256 seconds which is 
the size of the blocks analyzed. These features persist with 
the data block connection scheme and 144 point smoothing. Their 
origin is not understood. Some data blocks look like Figure 
3.19 with the end points having about the same value. However, 
some have end points as illustrated in Figure 3.20. The varying 
end point values are also shown in Figures 3.10a, 3.11a and 
3.12a. This may be a source of the large amount of hash between 
2 hertz and 10 hertz seen in the PSD's. 

After looking at the magnetic field data it was obvious 
a validation experiment was needed for the sensing system and 
software. 
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Figure 3.16. X Coil Magnetic Field, 17 August 82, 2241-2315, Local. 

Amplitude (nanoteslas : 1 unit/in) vs. Time (seconds : 500 units/in) 
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Figure 3.17. Y Coil Magnetic Field, 17 August 82, 2241-2315 Local. 

Amplitude (nanoteslas : 0.5 units/in) vs. Time (seconds : 500 units/in) 
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Figure 3.18. Z Coil Magnetic Field, 17 August 82, 2241-2315 Local. 

Amplitude (nanoteslas : 2 units/in) vs. Time (seconds : 500 units/in) 




Figure 3.19. 



One Frame of Data with End Points Having 
the Same Value 




Figure 3.20. One Frame of Data with End Points Having 
Different Values 



87 



IV. SYSTEM AND SOFTWARE VALIDATION 



A. EXPERIMENTAL APPARATUS 

The inconsistencies in phase between the time series 
voltage and magnetic field plots prompted an experiment to 
validate the sensing system and software. The sensing 
apparatus consisted of the X-coil and a Schoenstadt fluxgate 
magnetometer. The sensing axis of both sensors were horizontal 
and oriented North-South. Figure 4.1 shows the test set-up. 
Instead of an on axis arrangement, the set-up shown was used 
because of limited resources. The fluxgate magnetometer is 
closer to the source since it has less sensitivity at the 
frequency of interest than the coils. Recorded on the usual 
X-channel was the output from the X-coil. The Y-channel 
contained data from the coil after passing through an additional 
Krohn-Hite model 3321 lowpass filter with a 10 hertz cutoff. On 
the Z-channel was the output from the fluxgate magnetometer. 

All three channels were recorded simultaneously and a 64 
sample/second sample rate was used. An on-site chart recorder 
monitored the fluxgate and double filtered coil voltage output 
real time. The experiment produced two digital data tapes each 
approximately 46 minutes in length denoted GMTT1A and GMTT1B. 

Artificial magnetic fields were introduced into the region 
by applying sinusoidal currents from a Wavetek signal generator 
to a cylinder (diameter = 0.13m) wrapped with 150 turns of 
wire. Location of the source is noted in Figure 4.1. The 
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Figure 4.1. Test Arrangement 



applied sinusoidal currents varied in frequency from 0.005 
hertz to 10 hertz. For each frequency, two magnitude currents 
were applied to the source to demonstrate a linear response with 
amplitude . 

The fluxgate magnetometer output is correct in phase and 
amplitude to at least 1 hertz. The transfer function between 
voltage and magnetic field for the fluxgate magnetometer is a 
constant of 10 nanotes las/ volt . Thus the fluxgate provides an 
excellent check for the coil transfer functions, the RD data 
stripping subroutine and the digital fourier transform. 

B . EXPERIMENT RESULTS 

1 . Voltage 

From the real time chart record of voltage for the 
fluxgate and double filtered coil, we can check the computer 
generated voltages using VOLTR and VOLTS. Table 4.1 shows the 
chart record and computer generated voltages for the fluxgate 
and double filtered coil at various frequencies. Only the 
large amplitude oscillations are measured for the double 
filtered coil because of poor signal- to-noise , for the small 
amplitude oscillations. 

Table 4.1 shows good agreement between the chart record 
and computer output for the double filtered coil. However, 
the fluxgate magnetometer voltages undergo a constant 
transformation of 0.831 from chart record to computer plot. 

This may be explained by a problem in the data stripping 
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e filtered coil and fluxyate magnetometer voltayes, 
time and computer generated. 



algorithm. The present agreement between the chart record and 
computer generated voltages is shown in Table 4.2. 

Another more serious problem became obvious during the 
examination of the computer generated voltage plots. For 
approximately the first 20 minutes of data recording the 
fluxgate magnetometer was not connected to the PCM board. 
Figures 4. 2-4. 4 show the two respective sensing devices. 

Figures 4.2 and 4.3 both contain "real" signal while Figure 
4.4 does not contain "real" signal until the amplitude 
discontinuity indicating connection of the fluxgate 
magnetometer to the PCM board. Figures 4. 2-4. 4 indicate that 
"pick-up" is occurring probably between the channels at the PCM 
board. The amplitude of the cross-talk voltage is about 0.15 
volts peak-to-peak at the frequency of 0.1 hertz. The 
amplitude of the "real" signal on the double filtered coil is 
0.29 volts at 0.1 hertz. Thus 51 percent of the "real" signal 
is coupled to the disconnected channel. Coupling between 
channels of this magnitude raises serious questions about the 
computer generated data, especially for coherence. Discovering 
the amount of the channel pick-up if the channel is loaded 
requires further experimentation. 

Figures 4. 5-4. 7 show an interesting feature of the coil 
data versus fluxgate magnetometer data. On the channels 
containing the coil output we see noise on top of the signal. 
This noise is not present on the channel containing the 
fluxgate magnetometer data. Since all three channels passed 
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Figure 4.2. X-Coil Voltage, 26 April 83, 0931-0948 Local. 

Amplitude (volts : 1 unit/in) vs. Time (seconds : 200 units/in) 
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Hz Pick-up 0.1 Hz Signal 
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Figure 4.4. Fluxgate Voltage, 26 April 83, 0931-0948 Local. 

Amplitude (volts : 1 unit/in) vs. Time (seconds : 200 units/in) 



through the PCM encoder and decoder, the noise is probably in 
the coil amplifiers or the coils themselves. An experiment to 
investigate the true origin of the noise should be done. 

2 . Magnetic Field 

Given the dipole equation [Ref. 6] and the experimental 
arrangement in Figure 4.1, the ratio of magnitude of the 
magnetic field seen by the fluxgate to that seen by the coil 
can be calculated. From the geometry in Figure 4.1 the field 
at the coil would be a factor of 8 less than that seen by the 
fluxgate. This is because the source is a dipole and the field 
of a dipole falls off as l/R . Detailed consideration of the 
geometry gives a ratio of approximately 8.09. Table 4.3 
contains the magnetic field measured at the fluxgate 
magnetometer. Using the ratio calculated from the geometry, the 
expected magnetic field amplitude at the coil is listed in 
Table 4.3. 

A modified form of computer code LFVTC1 generated the 
magnetic field versus time plots for the test tapes. Table 4.3 
shows the computer generated peak-to-peak magnetic field 
values read from the computer plots. From Table 4.3 it appears 
that the coil magnetic field amplitude is in error by a factor 
of approximately three (with exception to the low frequencies). 
One can conclude that a problem exists in the software. 

Figures 4. 8-4. 9 show the magnetic field data for the 
second seventeen minutes of digital data tape GMTT1A. Figure 
4.9 is the magnetic field measured by the fluxgate magnetometer 
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Table 4.3. Measured and Computer Generated Magnetic Field Magnitudes 
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Figure 4.5. X-Coil Voltage, 26 April 83, 0948-1005 Local. 

Amplitude (volts : 2 units/in) vs. Time (seconds : 200 units/in) 
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Amplitude (volts : 1 unit/in) vs. Time (seconds : 200 units/in) 
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Figure 4.7. Fluxgate Voltage, 26 April 83, 0948-1005 Local. 

Amplitude (volts : 0.5 units/in) vs. Time (seconds : 200 units/in) 



and Figure 4.8 is the magnetic field seen by the coil. For the 
coil we see discontinuities every 256 seconds corresponding to 
the block size of data. This indicates that the problems in the 
data block end points, as mentioned in Section III, do exist. 

It is not clear how to deal with this. 

The first 300 seconds of Figures 4. 8-4. 9 is expanded in 
Figures 4.11-4.12. With the expanded data a phase difference 
between coil and fluxgate can be determined. By checking the 
phase differece, we can determine whether the phase independent 
transfer function is introducing error. For the 0.1 hertz field 
plotted in Figures 4.10-4.11 the coil lags the fluxgate by 
approximately 72 degrees. In Figures 4.12-4.13 we see that at 
0.01 hertz the coil sensitivity has fallen and a reliable phase 
difference between the fluxgate and coil cannot be made. 

Figure 4.14 shows the relationship between phase angle and 
frequency . 

The above experiment shows the advantages of having a 
real time "ground truth" record. With this, one can avoid 
"working in the blind". Computer products will have a check. 
Thus, a multi-channel chart recorder should be purchased to 
graph in real time the x, y, and z inputs to the PCM board. In 
addition, electronic circuitry should be developed and 
incorporated so that the chart records represent at least a low 
pass filtering of the coil voltages. Ideally, conversion of 
voltage to magnetic field could occur in real time at this point 
given the appropriate electronics. 
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200 units/in) 
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Figure 4.9. Fluxgate Magnetic Field, 26 April 83, 0948-1005 
Amplitude (nanoteslas : 10 units/in) vs. Time 
(seconds : 200 units/in) 
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ure 4.12. X-Coil Magnetic Field, 26 April 83, 1023-1033 Local. 

Amplitude (nanoteslas : 0.5 units/in) vs. Time (seconds 
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Figure 4.14. Phase Lag of Coil with Respect to Fluxgate 



V. MICROPULSATIONS 



A. THEORY 

Currently there is no theory that completely models the 
observed micropulsations and geomagnetic noise. It is 
generally accepted that the generating mechanism for 
micropulsations is an interaction between the solar wind and 
the earth’s magnetosphere [Ref. 5]. The observed 
micropulsations are classified by their period and regularity 
[Ref. 5]. Table 5.1 shows the micropulsation classification. 

The smoothing algorithm applied to the magnetic field data is 
designed to highlight micropulsations with periods between 10-45 
seconds. Specifically we are interested in micropulsations 
classified as Pc 3 . Pc 3 micropulsations are generally regular 
and have a more or less distinct period. Thus in the PSD plots 
we will be able to identify the micropulsation by an apparent 
anomaly at its frequency. Micropulsations have amplitudes of 
approximately 0.5 nT . This small amplitude is very small 
compared to a main magnetic field value of 50,000 nT , hence the 
term "micro" pulsations. Geomagnetic fluctuations are a 
combination of micropulsations and a general background noise. 
Micropulsation events grow out of the ever present background. 
The time series magnetic field plots will provide a means of 
obtaining the frequency of the micropulsation. 
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Table 5.1. Micropulsation Classification 



Notation 


Period (seconds) 


Pci 


0.2-5 


Pc2 


5-10 


Pc 3 


10-45 


Pc4 


45-150 


Pc 5 


150-600 
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B. MICROPULSATION DATA ANALYSIS 



The first 1200 seconds of data after a 60 second tape 
advance on digital data tape GMDT 11, 17 August 82, 1301-1408 

Local, may contain a micropulsation event. Figures 3.10-3.12 
show the first 34 minutes of GMDT 11 after a 60 second advance. 
To determine whether the regular pulsations we see in the first 
third of Figures 3.10-3.12 (disregarding the suspicious "cusps") 
is a micropulsation, the scale must be expanded. Figures 5.1- 
5.3 show the regular pulsations on an expanded scale. From 
Figures 5.1-5. 3 the period of the oscillation is calculated to 
be 17.2 seconds. Thus this section of data contains a type Pc 3 
micropulsation. 

Executing the computer code LANDXYZ [Ref. 1] on this 
section of data, we generate the PSD, coherence, degree of 
polarization, and ellipticity. Table 5.2 contains the 
amplitude, frequency, coherence, etc. for the micropulsation. 

At the frequency of the micropulsation one would expect a 
good coherence, degree of polarization and ellipticity since 
the coil system is oriented arbitrarily with each coil sensing 
the micropulsation. The observed micropulsation exhibits a 
high coherence as shown in Table 5.2. The observed micropulsa- 
tion has a high degree of polarization. Micropulsations are 
generally considered to be left-hand elliptically polarized 
[Ref. 5]. The observed micropulsation is elliptically 
polarized . 

If the errors in phase and amplitude that were shown to 
exist in Section IV are constant, then that should not be a 
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Figure 5.1. X-Coil Magnetic Field, 17 August 82, 1305-1310 Local. 

Amplitude (nanoteslas : 1 unit/in) vs. Time (seconds : 50 units/in) 




flOO £00 £00 100 000 
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Figure 5.2. Y Coil Magnetic Field, 17 August 82, 1305-1310 Local. 

Amplitude (nanoteslas : 1 unit/in) vs. Time (seconds : 50 units/in) 
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Figure 5.3. Z Coil Magnetic Field, 17 August 82, 1305-1310 Local. 

Amplitude (nanoteslas : 0.5 units/in) vs. Time (seconds : 50 units/in) 



Table 5.2. 



Observed Micropulsation Properties 



Frequency (Hz): 0.058 

Peak-to-Peak Amplitude (nT): Amplitude ='^X 2 + Y 2 + Z 2 = 0.37 

PSD (dB) (Ref lnT 2 /Hz) : X-Coil: -3dB 

Y-Coil: -4dB 

Z-Coil: 1 . 5dB 

Coherence: X-Y Coils: 0.99 

Y-Z Coils: 0.99 

Z-X Coils: 0.99 



Degree of Polarization: X-Y 


Plane : 


0. 99 




Y-Z 


Plane: 


0.99 




Z-X 


Plane : 


0. 99 




Ellipticity: X-Y 


Plane: 


0.45 


(left hand) 


Y-Z 


Plane: 


0. 10 


(right hand) 


Z-X 


Plane: 


0.30 


(left hand) 
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serious problem since we are dealing with a single frequency. 
Thus , the observed micropulsation exhibited definite structure 
decidedly different than the background. 
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VI. CONCLUSIONS 



A type PC-3 micropulsation event of peak-to-peak amplitude 
(0.37 nT ) and period (17.2 seconds) was identified in a section 
of data. At the frequency of the micropulsation, the coherence 
between the three orthogonal coil sensors was about 0.99 for 
the xy-pair, 0.99 for the xz-pair, and 0.99 for the yz-pair. 

The degree of polarization was 0.99. The ellipticity was about 
0. 30. 

The origin of the regular "hump"-like structure in the 
PSD’s of previous coil data was found to be artificial and not 
representative in any way of natural phenomena in the 
geomagnetic field. Previous conclusions based on a frequency 
domain analysis of coil data characterized by a regular "hump' 
structure (strong or weak) in the PSD plots are known to be 
unsubstantiated . 

Compared to a PSD frequency domain analysis of the data 
set, an analysis of geomagnetic field fluctuations in the time 
domain easily identifies data sets containing a micropulsation 
event. A double running average routine developed and applied 
to the time series data sets was found to be helpful in 
highlighting the micropulsation. 

Detailed conclusions about the geomagnetic field from the 
coil measurements require a balance of analysis in the 
frequency and time domain. 



118 



An accurate reproduction of a time series plot for the 
geomagnetic field oscillations will require modification to the 
existing software. 

Comparisons in the time or frequency domain between any 
pair of coils in the orthogonal system might be misleading as a 
result of cross talk between electrical inputs. The suspicion 
of cross talk is based on the observed pick-up on an unloaded 
input of the contents of an adjacent loaded and active input. 

Further micropulsation studies should continue. This 
initial investigation of the coherence, polarization, and 
ellipticity features of micropulsation events shows promise 
that with more data, one might suggest novel physical 
mechanisms and models for the origin of micropulsations. 
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APPENDIX A 



VOLTR COMPUTER PROGRAM 

Appendix A contains the VOLTR computer program that 
provides unsmoothed time series voltage plots. VOLTR uses the 
subroutine RD to strip voltage data off digital tapes. The 
integer voltage data is normalized to + 5V and plotted using 
the DRAWP subroutine. 
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APPENDIX B 



VOLTS COMPUTER PROGRAM 

Appendix B contains the VOLTS computer program that 
provides smoothed time series voltage plots. VOLTS uses the 
subroutine RD to strip voltage data off digital tapes. The 
integer voltage data is normalized to + 5V. A double 144 point 
running average is applied to the data. The DRAWP subroutine 
plots the data. 
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APPENDIX C 



LFVTC1 COMPUTER PROGRAM 

Appendix C contains the LFVTCi computer program that 
provides smoothed time series magnetic field plots. LFVTCI 
uses the subroutine RD to strip voltage data off digital tapes. 
The voltage data is normalised to +_ 5V. The voltage data is 
fast fourier transformed using the FOURT subroutine. The 
transfer function is applied and forward transform taken. The 
resulting magnetic field data undergoes a double 144 point 
running average and is plotted using the DRAWP subroutine. 



131 



OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

^cMmM"mor^-coa>Or-jrMm<MnOr*-coa '0 — 4c\irnM"tn or^coO'o— 4<NcO'4-^'Oh-coo , '0— 4 cmcOvJ- in or^-ao 
o O O o O O O O o — * — *■ — »< — < —4 — < — < cm CM <\) CM cm <\ j eg cm cm cm m ro rn m rn m <n ro m cn ^ M* >r m* ^ vf <r m- 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
<<<<<<<<<<* 









<LU 




o 






Q 


aco 




2 






2 


— ^ o 




— » 






< 


QUh 












LUIU>— » 




< 






< 


2Q.Z 




ac 






1— • 


a: loo 




LU 






< 2 — 


o < 




2 






QO(N 


0 . 0 : 3 : 




LU 




X 


— O' 


OOLU 




CD 




CO 


H- 












O<co 


<02 




2 




0 


CLX~— 


a: cllu 




•—4 


0 


CM 


20CCMCM 


^ LO 






II 


11 


— QOO 


OJJ 




a 


00 


a 


LLOl— * 


LLh-CX 




UJ 


00 


0 


JJ — OC2LLCO — 


u- a. 




OO 


< 


• 


q-mlu— *~vO 


wOLU 


-r 


O 


— 1 


2 


<o^q 


LUQC 


jj 




0 


0 


4 c£ JO H- I_n 


2h 


> 


LU 


•* 


♦— • 


co 0 lu co fvi n 


—a: ^4 


< 


a: 




0 


x*-^ >— »m *-sO 


<UJ^- 


•* \ 


< 


0 


'jj 


QU.Qa' O — w 


SG>rv| 


CO 




x 


-1 C£ 


Q£ h— H» JJ —4 C\l *— 4 


02 


UJ H- 


m 


0 


< - 


a. •* oo^*o> 


OGDO 


> a 


LU 


(M 


0 - 


—0 x:— 4 .^ 


02 


< H- 






Qll 


<nj JJLUOICOM 


> < 


•» W 


h- 


U 


-JUJ 


H-c^oosiO'— - 


— o> ► 


PsJ 




^_ 


il a: 


<£ 0 *— * 2*—**— 


>0 2-1-4 


JJ • 


r— 


00 


H" X 


IJOO 1 — -> OcOULIJ> 


>1- M 






00 - 


W UJ — rvjrninOC^M 


< H- 1 






JJ Ql 


ONt 1 CM -nmoo •* 


-00 X 


o • 




— o< 




—42 o 


2?- 


00 in *»2 


— 1 •* UJ — 4f\| sO 


— a:2> 


'JJO 


— <o 


^ . 


ooc - 


.0 — 00 O C0O N — ' 


CMLLLUrvl 


>o 


— ••CL 


jj' 




< CM >- < ^— 4 — 4 


LU "> *• 


< - * 


n- !— 


> II Ooo 


JJO<1—OCO>- 


2^-’JJ— * 






•ajoo'U— • 


j^-4 0^-JLU— ^ 


— 2 X 


O 1 


0<0— x 


1— UJ>0 


CO rc O CC -4 1 M H- LU LU PM 


Zlos *X 


h-> -J • • • 


</)2 ,y >- 


z^<>u.> - 


•*00 QC 


GZOOO 


a:<aioo\so 


• 1 


—>• 11 


•*-g— - 


— LU< • 


uoo-^ - 


•*• UJ — •* • •& 




-JOOQ 


>-coh-— z: 


-orgco - 


— -f— • *OOM* 


— * 


» • JJ 


0 H Zm^cocoq-QO 


fv| •»»^ CM • 


-4- OO-il* ^ o 


OCL>— ^ 


lu— x lu ocMinm lu 21 ^ 




— CQQw_<fM-m 



— «cm~J • 
000.2 
— < II GC 
O *LU< 

cm— * 20 - 

Oi< •**• 

-<> 2 : 0 . 



oocMuuoo^-4CMnmc*c<*— 
-OO'-juuco— 4«ogo < 
0 --o.tX'— cO'-'—oocs: 

--400 00 ZD-JJ^— 4— 4>-|— q£ 

^~QtiiS^Xh«0 
2 XuXhXNNqJQii 

X rsirsiivifX 2 

LU >• <LUh-« 



OOh^ I LU|— _J-fr # OOOOsO 
•■oO 00 cm GO X — -J QC cl vf rn m —4 
in -z— ' cmcjcmi— - coco nox 

*— «>»0— — r—4 — » ••omen — <— 40 
-XOcocax— o- .jj so s 0 'n* N s uli 

<y ^ I— ' LLJ CO h— —“4 —“4 —4 1—4 QC 

jLUH < tNN>h*LL 
— Oh •O 0 >»U_I s nIM » 

xxoo>t o»— 2 - - -lu 

22 tt <h-uj _ac 



^OQOQN«00>-JZZZ2 Q 

HG^Q> - ^.OJOQQO- 3> - — 

02 ^ a: x ( o<* , ^'~ 4 HsHH NHt— </) q^"oCoc: _j j^ulixnx>^ 

Q II (XXXUUVO^t^oOoOoOOO*. — ‘lUuJoOO'JJM^COO-XXMXX*-- 

ooah-cxj u-zzzz 2 looo CD# # x 

OC C£ f /}UJC£CLiJJ»— OyjLUUJLUC02— IUJOO , JJ_J_J__J’"" , rY'-U<<<<< ^-»-4(/)tO00^2X>N^-4(\jrOvf 

-■ u_ <r:nci~c ct~<<<z>a: iTi—^ O 11 H 22222 XX Z’-uiulj'Ulu 



0000 

• • • • 

000000000000 

II II II II 

X>NhOOOOOOOO 
II II 11 II II II *1 II 



;zzszzi<ai<H«h 
< ooz<oi— *-•— -♦■——i- xo> 2 ujuj 2 ujujluo<>-<<<< < n 'T^oaauoo^J>>» 
r^ZZo • O OQG O — 4QCo: — QCc^CC'-U QQQQQX^OOUU00t0^00<<« 

I— — OC'JJH- 
><1QXIX 
izxaja 
o-* * a. 



o 00 



ouoo 



00 



132 



OOOOOOOOOOOOOOO OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
o cm cn m vO > oo O' o c\j m ^ in o h* oo as o cm m >4- in h- co o o »-* <\j m >*- vO > 00 o o — cm <n ^ ^ -o 
^ininLnininininintnin'OvOO OsC 3 ^vOsOvo>>r^>>>>>f s, --H-oococococococooococoo ,> oC v ^oO'o-' 

0 000000 00 00 00 000 000 0000000 000 0000 000 000000 0000 00 
00000 0000 00 0000 0000000000 00000000 000000000000000 
000000000000000000000000000000000000000000000000 



<<<<<<<<<* 



C<<<<<<<<<<< 

: 2: 2: 2: 2: 2: 2: 2: 2: 2:2: 2: 2: 2: 2: 



X 

lu 

H- 

UJ 

V 

a 

H- 

UJ 

2 

o 

< 

2: 

LU 

X 



X 

< 

> 



OO 

o< 

ox 

s 

— 'a: 
0.0 
O 



h-O 

— * # ■— I 

^ t ^ 

• r\| • Q< 
O-iO X 

NSO I 

a 

i^aa\ *o 
3 )OQh ouj 
— '^^Gl OQ 

mhH-o^ 
ao— ’'—xa — uj 
CMoOlOH-^— 00X 
• OO# oOOh- 
000^00 
II II II H-O || OUJ 
— OO • II 
Ql O m JD O O — 

iau 11 ua 

I— OOQOO 



rO 

in 

in 



UJ 

x 



o 

X 

t— « 

H- 

x 

< 

1— 

0 

X 



> 

< 

0 

JJ 

a 

<jj 



o 

< 

UJ 

> 

X 

JJ 

00 

o 

UJ 



JJOO 

— 00 

LL> 



X 



o 

IJJ 

or 



o 

0 



00000 • 

•* • • • • o 

-<0000 II 
II II II II II — w< 





0 


UJ 








• LU 


2: 






a 


QZ 


< 




> 


H-X 


LUX 


X 




0 




o< 


LL — 1 




XX 




<x 


1— 1 




XUJ 


<00 


xo 


za 




>03 


or 00 


LU 


00 




OLLGr 


OUJ 


>x 


Ql 




X LU 


00 


<0 


> 




UJUJ X 


oca 


< 


LUO 




ooll 


aa: 


UJUJ 


X 




o< 


a 


X 


H-< 




JJI— LU 


UJ 


X 






XJI 

U,QH 


XUJ 


□0 


x< 




H-I 


h- LI- 


a a 




> 


H- 




x 




u a 


00 


OO 00 


LL UJ 




ZOH 


LUO 


Uh- 






< 1 — 


-IX 


ox 


<H- 




2 


CD 1 -* 


X*- 


H- 




JJCX 


<h- 


luq 


<H- 




SHQ 


x< 


3 CL 


OO 




•— < LL 


LULU 


3 


a 






CL 


UJ< 


UJ 




X 


ouj 


OO h- 


XLO 




'XI 2 < 


>x 


< 


1 — ►— 




J-OX 




<3 






HO.H 


h- > 




'y)Q 






XO 


<GsJ 


3 <J"> 




yO|— X 


'U 




< 




JGLU 


x< 


—4 


OO 




H-X — • 


JJH- 


LL 30 


XX 




<xx 


>-< 


a 


< 




x-o 


<0 


00 


0 




-u 3 


t— 


X — 1 


JD •» 


X 


X JJ LL 


oOLL 


UJ< 


u 


X 


ill 


0 


TJO 


X G 


— • 


5 HH 


T* 


zcr . 


H-O 


•» 


00 


H- h«B 


OUJ< 


•^•O 


O 


XO< 


— X 


X H- 


3 


LU 


CUJLL 


^0 


>< 


O 


X 


•— •f'Nl 



a 

ox 

x< 

OUJ 

xo 
lux 
— x < 
azau 
im x *a 
Ch -*~*o< 



LU_JQ 

Ih 



02 

x> 



h-~x 

oao 



c\j •» •* 

»— <0 

XX # II CNJ —I > II 

LULU< UOwU'O <t —< 

^2Z2Z m OZ OU^GZ) II O^lOoOoOcCO'O 
<^Ooo xXlu -J Qlux am 

^,w^.wrg^uj|- i| i—in mJ LUO Oo- UJOvj 
r0^^^^UJh-X<O II m_ju-<00 II MllOaaCLO 

ZZHQOJJ jza 11 X I'UZ X 

Q^^rNlM— a lo>— Q<OLLX ZOh O.CQ*- — *Q 

ONNNNhO ^ M QOO^XU,Q 



l-XLL 


t— < 00 


O 


LU<L 




UJO 


0 


LUO 






oox 


XUJ 


X 3 


XX 




f— XOO 


UJX 


< — « 


oao 




xoa 


1— 


X - 


XX — 




UJOX 


X 




— V— 




00 0 


U(— 




^a< 




UJLUO 


uo 


- v^- 


axx 




XOUJ 


00 


1-4 0 ' w ^— 


_j<< 


X 


azoo 


-J 


♦! CMXXX 


0 x • 


— » 


UJUJ 


O OOLU 0 ^ ' 


Luaoo-jjx 


H- 



Qazoa 11 u 11 oll>x-- ■ < 

— X X~~ * — x <Q.<(\JO 

ujx<o 0 - 00 — lux xo_j 

TH-TOjO-JOhl^ZC-iU. 
h-OOO _J — — "■'X)— <— QO II 

a<x>No nx 

QOX>.MO XLL 



co 



in 

n 



o 

o 



00 



00 



0000000000 



0000 



133 



OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 3 
go o —i f\i co m vo I s * oo o o — i cm co m 3 I s - oo o o cm on m so I s - o o ^ <n co <j- in 0 I s * 03 o o — < <m co 
o 0^ o o o o o o o o o o •-» -4 ^ -* -* cm cm cnj r\j eg og c\j cnj c\j cm m co m m on m on cn m ro v$- st > 4 - ^ 

O O O >” 1 ^ r “^ r-4 F— 4 p— 4 p“4r— I *— 4f-H »— 4 r— 4 *“H r— 4 P-^ *— < P- ^ p— 4 p-4 f4f4f4f4f 4 F-Jf>4i->4r-4rH t—i p-J r— I fW —4 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

oooooooooooooooooooooooooooooooooooooooooooooooo 



UJ — 






H-» 


IN 






COLLI < 


1— fNl 




a 


<H-X 


— 




h— 


X<fa: 


LL 




> 


zao 


OH 




UJO 


<3LL 


X 


h- 


OX 


oox 


XUJ 


3 


<UJ 




ax 


0- 


^-3 


< 


»>»' 4 Q 


X 


03 


COXLU 


<1—0- 


—a: 


Qai 


<*-HLO 


H05! 


o 


>QC 


a < 


<ujq 


lull 


wtL 


oozi 


030 


X 




LU<Q- 


a 


1— - 


IQ 


• H— 


joa 


H- 


UJ.UCMGLOLL 


•— »a. c 


XQC 


HZ 


ZQh. 


a o 


03 


gocx:h-<o 



OUJHH 

t th 
>h-q: • 
luqc: 
UJOO >0 
Im 1— 
I— uu 
• zn ~u 

LOLL H- L> 
—•H— 

- OQ 



a 

KUL 

LL - 
LL 

x5 

< 

a. 

X3 
x a 
31— 



> 02 ZUJ 2 
COLLUJ mQ 

cos:a>— 

UJXUJH- H- 
X<1 — O0U 
H-cC«n2X 
H— ♦— ■ *— O 3 
> lO— I — 0 
-Juj Of» 

n.Z^>uoz 

CL r— < *ZUJ 

< V)3LL 











>-0<LU 




OtLcC QOH-lLO 


















^X — 




0J3 LOI— goc X 




33 


— — 






• 


• 


• 


• < — LL 




CL — — H- X'UQC< 






r^os— 






O 


cn 


o> 


X 




jj X^^ZZUJ>iua 




OLLt-L 


O>O N U0 






o 


CO 


m 


~ *xo 




oO X C^LLIO XLL^— 




->!• 


• • co 






o> 


3 


Lfi 


<<r^— • 




H-h- 300 X— »XOlo 




U-OvH 


oo • 






f-4 


CM 


M 


h— H* H- 




X 332UJHOOZH 




* 3 3 


corn m 






— — 






<<H*UJ 


— 


lulu *• •* -h-oo <x 




sJ-lAM* 


1 1 cO 






3-5 • 


• 


• 


ooxx 


— 


XUJ *000<ZJZ^UJ 




f-4 • • 


33 | 




— 


— ^Lf> 


l n 


in 


LUO 


3 


LUO02 •* •* ►h-*3CQOl— QC 


f-4 


(M •r s -r— 


cncx.cc.C3 


<r 


H- 


h— • H” 


* 


# 


aox< 


— 


h— 3 H f— 4 >— 4 1/) LL JJ 




cO | I 


U-LL CC 




< 


<< — 


— 


— 




O 


< •— 1 I 1 | OOI — OOLL 


a 


a i 


O* * LL 


o 




03 • 


• 


• 


ooq-o 


UJ 


i — <z h- •* •* »*ll a:— *o — *ll 


H- 


i— mcocM 


H-3 3 


H- 


a 


JJO 


— cn 


— O 


^UOZHI 


GC 


'vOh-<^ >— * ■— 1 0 O IZX M 




• • • 


aoo o> 




jj 


LL LL O 


— CO 


—CO 


— 1 1 GO 


LL 


<x: - - ILt— 3h- O 


a 


3in~ir- 


OOs^nJ- 


3 


a 


> * 0>303lA-5X>UO 


— 


GC3QCXXXXO LL 


o • 


OOCOH 


o • • • 


O 




h— LL — < 


—CM 


^ryj 


— -J 


O 


3 a - - •‘OZ J * < 


— CO 


— F— 4f— 4 f— 4 


— 


— 


•X 


<< 1 


X 1 


> 1 


mluluX< 




Oi/)lUX>NQ< «o^* 


•<M 




• W >MI 


in 


CM U- 


f— t— — » 


X — 


>• — 


niihh 


O 


LLuJ2X>NUq:<UJ^ uj 




i^INNS 




• 



co— X- 



II 



5— 3 — “5 — 1— H-3'v 



2ZU 






UJ 



MMww^CQf-h-LL^Jt/) 

«v. .*0 0—0—0—0 OH- I— a: H-H-H- XO »0 »000 •— IOO 

•f4Qcx<><n<^oO(/i f-4< x o a: ac cx: qC h- —* Q^h uj— o— -- o— — — 

— i l| l| II XUJ>UJ!V|UJ ►-■»—« | LUUJ II II LUUJG0UJ333X3 <2 *^OjX JX>N _IX>N O 

,l -5 .-wr£_a: — a: II3 ^3Z XQOGLUjZoia: «<3 n uj *xco «x>-mco «x>-mco • 

r— LL 33 l| II II II II 1 | - • h-H- 2 3 ^ LUl— LLLLlL 2 UJ — *H-<OUJOaC 3 II Oil II II OH II II 3 

<<0-— — — .— — XMC£ h-.p-4^-p— o SlCC •— '< 3^0 LLQ£ — OC£ — — <3cC~*~ — Oa£ 

H-l— (MOnT33— 5333 XIVJOX!— <NcMH-X^ 30000LL:£0-* <LUO ll LL — I h— U I I — 1 1 — U lOOl— LL 

OO 2UJ'-— • ZOZ OZHHU U-JI | dl^IUJ 3 — ~ — 

ujajUH-.oixx>->rsjrNsi ouixa <<<i— ajQh--* a.z'jau.x jLX>Naxx>NUiL 

QQQHU.XX»NN UQXU OOO QO.^XO'-.X>NlO > -'X>No^ 



o 

(M 



CM 



ro 



cm 



oooo 



ooa 



O OOO OOO 



134 



OOOOOOOOOO OO 30000000000000000000000000000 0000000 

m o co ct' o ^ <\j ro in h* oo o ^ cm no co on o ^ cm ro o ^ o o o <m ro sf * n o co o o ^ f\» 
<f <• <f < <f in in in m in >n in in in m uo vO o o 'O o * o 'O r** i^- I s - js. r^- h- I s * r*» oo oo co co co co co co oo oo o> o> o 

4 *— -4 »— 4 r— 4 t_4 «— 4 rHn4 *— «— **—4 *-*«l r— 4 .*«4 •— r— 4 r— i«4 r— 4 r— 4r— 4 *"“4 r—i r— 4 p—4 f— 4 r— 4 r— 4 r— 4 1—4 «— 4 •‘■‘4 —4 ,—4 H»^f-4H r— 4 «— »-^ »— 4 rHr-4 

000000000000000000000000000030000000000000000000 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO'OOOOOOOOOOOOOOOOO 

<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 



o 

ro 

co 

O 

o 

ft 



• 

00 

01 

2: 

o 

o 



rvJ 

rvj 



o 

o 



o 

o 

oO 

a 

o 

ft 




o 

LO 

a 

o 

ft 



X>*NIL 
X>NH 



< 

o 









— ■* 




O 




■S^ 


r- 


<— 4 




-) 




>0 


^0 


0 




at: 




NO 


■O 


</) 




at: 


— <n— 




vf 


O 




•0 


— 00 


— Lfivt 


— 


O 


— — — 


0 


ro *C0 


• in %r 


• 


ft 






• O * 


0 • • 


O 




oC&cZcC 


LU 


0 -- 4 N 




f 




aooa 


O 


1 1 I 


O | 1 


O— 


— • 




O 


000 


aroo 


a: or— — 


> 




O 


ai^cx 


LL CCCC 


LLocar O 


> 


— — — r — 4 





co co m co 



X>rviU. 



2COOOcO<^0 

►ajaoaJcQ 
— <<<< 
|| ouuu 



u. 

OvO 
CO CM 

OO ONX>N^- 

**oOcOoOoO 
JCCO + + 4* + 

sT ococooooo 

» 00 t— II— I >-M— I 

O II 



4--Tvf vtO' 



U-J — — x O^aToro: — — — - II 0000 O II — w — OOOO l| 

rM x UJXUJOOOO II X>rvlLULUrO il II II II JJL- *oO«^— • — • — iUJcOlOoO^oO 

♦ X>NJ0DXSD^DUOQO-)X>NhDM^ OX'-U^*X»KD II !| II II — 

II II 11 or II II II o 11 II (| (I 2 1I2 U.lLU.LL 11 II II M 2 c0rcunr02UJ2 rv|rx|rMrvl2X>fM^ 

*0^0^-*^ *-*vO— — — — •-•2 *o u u u »| — i— ^v— f—h» 

i _j IU. M 1 |_i^— LI IJ Jh JhL _J|— . — 1 j_ Iin~3 m3L*iriw^ *-» — . j— ^v(5X>NLH cOcOcOI/Vq 

— w—™*-. 2 — 2.J-J-J-J — — 2 — — — — 2'UJ — SIIt22Z22 

X>NQa.X>NUlLX>NQXaUU.U««CX>NliUaX»HUI*ia^033UUGC f J'J 

X>NO^X>NO^X>nIOXOUHOUUUOQX>NLOQNNMNUHmO^V1^00OU(JOJG 



in 



ao O' 



in 



in 



vO 

o 



135 



OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO oooo 

m st in so r*- go O' o — i cm m < in o I s - co O' o <"\i co <t- in o r** go on o ^ c\j m > 4 " m o go o> o — « cm m vj- m o r* go o' o 
O' onqn O' O' o> O' oo o o o o o o o o ^ r\ic\i c\j cm cm rg <nc\j rsj rum rnrnm mm mm m m <■ 

^ ^ — « ,-i cm cm cm c\i cm cm eg eg r\j cm cm cm eg cm cm <M cm cm cm cm cm cm cm cm cm cm cm cm <M cm cm <m cm cm cm cm <N cm cm cm cm 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO'JOOOOOOOOOOO 
<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 



a 

a 



O 

2 



ID 

cc 

LLi 



^rgr^vt 
'Ll '-U LULU 

>>>> 
<< << 
I I I I 
X>Nh 



< 

o 

3 

LL 

■fr 

o 

» 

'S\ 

eg 



CD 

3 

a 

a 



uo 

x 

2 

a 

LL 

'X 

UJ 

Q_ 



a 

ID 

o 







cOcOcOCO 


< 


LL 










2222 


3 


a 










aaaa 


3 


• 










0000 


LL 


C 0 < 




x> 








* 


UJ!— 




XX 




X>Nh 


4 . + + + 


1 — 


2 < 




33 




IIII 


— 


< 


~Q 


00 


1 /D 0 O 


LOOOcOi/) 


33 X 3 


C\l cOcOcOCO 


t\JK 


3 




+ + 


»— <• i— * *— < *— < 


vJ-LOlOOOLO 


• • 4 • O' — < H - 4 1-4 »-H 


0'3 


ai 


m 


— ■ •— * 




g- + + + + 




— <UJ 


OX 


in 


sT 33 




rH —• 


g-g-sT sf go *— ••— < 1 


GOO 


Zh cm 


3 


-r + + 


X^^H* 


••CO 1 /) CO GO 


~X>->|- 


•*>— 


i—i » 


•» 


-*^00 


rgiNgNJhg 


*““• M 


SS\S-^MNNN 


-I II 


2 Z-i 


1 -^ 




II II II II 


11 


X>Nk II II II II II 


II — 


30 H 


11 


— J r — 1 




UJ I s - LU O O O O 00 ^ — « — 1 — l UJ s x DC 2 : GO «— — * — — 


LLJ LU min LU 


33 CM 


cOOOOO II x> 


sj-sTvts^^xmz) • • • *^x>>>— vj-< g- — * 03~-«-« *3 


33 LU 3 


1 — « 


• • * • 3 >X M 




X ZOOOO NNNrJZ^ oOcOcO n-« 


+ 22 — +2 


200 


OOOO •'MnJ 



si- ii ii ii ii in H ii 

r-x>* 3 H- r*-x> 
o 2 : 3 : 312 : 3 : 2 : 

»>aNNNN<ruuu^^u 11 ox>u 11 aDJJJcn 

UOOLOLO(/)»/)Qi/)CO(/)l/) 0 < ?<<<CDNNNNt-HOOQh' , * , U^OH<QOQLOlOiOl/)OCOl^ 



11 || || 11 00 11 11 ii ii 11 11 11 11 ow—~ tc\im»-»^ t - , LL<^n 

— i— H-H- X>NI- sOX>*3H-h-~<CMm<J-'0^^-^^*~*l*~H-0'LU>'-»l— ’ + i~ xr- 
X>>^- II 2 ZZISS XXX‘X2LUUJLUUJ X>>H- II zz 3: 11 222UJUJ 

rsj 

“M rg 'M r*g >— « 



{S. nQ 

sO m 



GO 

o 



0'f s - 

om 



O' 



o 

I s - 



JU 



136 



^cNmv$-mof^cooOr-HCMco>3-m o o> o*-< cm ro in x r** co <?* o-^ <Nm'tinof s -coc>o^ cm in 'Or^o 

<f N^vt^n in mm min irun la m soooo^oo oooor^-r^-r^p^r^r^.(^>-^-r^ ucoqocoooxcoooco 

cm cm c\i cm c\i cm c\j c\i cm cm cm cm cm cm cm cm cm cm cm <m <nj cm r\l cm cm cm cm cm cm cm cm cm cm cm cm cm cm cm cm cm cm cm cm 04 cm cm cm cm 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOO'^OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 

xxsrxxxxxxx 



X 

o 

LL 

CL 

LU 

O 

CL 

a 



Q- 

O 

I 

JJ 



CL 

IS 

LU 

X 



3: 

U 



a 

LU 

X 

f— 

o 

o 

X 

LO 



> 

CL 

< 



O 

u 

CL 

ui • 
uiQcn 
IULUUI 
Oh-D 
alhj • 

zo<* 

-JI>CL 

ooa. 3: 









f— • < 


—ft 


—ft 


— ft 




— 


—ft* 






O 


xujcdx 


0 


0 


0 


cn 


O 


n 






U 


— 'Oh- O 


i— 


h— 


(— 


f— 


f— 


}— 






u 


0 X- 


X 


X 


X 


X 


X 


X 






— 1 


0. 3 - 


•* 


•ft 


— 


•» 


•ft 


•» 






LL 


f— LU 


CD 


U 


X 


u 


X 


u 








LL OX 


f— 


f— 


V— 


f— 


I— 


f— 






O 


OUJXO 


— « 


— 


— 


— 


— 


— 






— 


o<o 


•» 


•ft 


•» 




•ft 


•ft 






h>» 


cL ^ U 




ft— I 






1 — ^ 


i —4 






'U 


u<~ u 


X 


> 




’r— 


X 


> 








0 x O 0 


-M 


'vl 






•vj 


-nJ 






O 


X l-X 


fvl 


rM 


•vl 




M 








< 


^)O 0 i— lOL 


•» 


•» 


•» 


•» 


•» 


•» 






X 


• 2:0- 


CNJ 


CNJ 


CM 


CNJ 


"NJ 


CM 








— ( UJ UJ 


U 


U 


LU 


UJ 


LU 


UJ 






1 


+ OOCOOZ 




X 


X 


X 


X 


X 


SIS 


• • • * 




UJ 


►*>«« 


I-— 


— 


1 


— 


— 


OO 




CQ 


V— 2 — •!— 


UJf— 


'JUf- 


Uh 


Jh 


LUf— 


LUh- LU 


0000 


^ vj- 




<'-*— »~2 0 


_J - 


u - 


u - 


— 1 •* 


U - 


-J - — J 


+ + 




LL 


1 — X O 




b-00 


f— lO 


I — 00 


h-v >0 


(S) f— 




"n, n ^ s v > s 


O 


uxoux 












— 1 — — — 


00 


X>rvlt— 




uj u CNJ u 0 


l- a. cnj i— a. m »— a. 


nT i— CL 


*— < <f—Q.CNJf— CLrO^™ 


4 * + 


Z 2 .ZZ 


f— 


Oh- 00 





— 


— X 


— -'X 


——2:— 




OO 


OOOO 


O 


NUJQU-^ 


00— » 


UO— ujO — 


xo— 


UO— cuo— cuO 




UHSHSHS) 


U 


• Of— 


00 < 0 Q-< 0 cl< 0 q. 


<OCL 


<OQ.<Oq. <0 



— < I! II II II 
^> 1 — LU'"'— JJLL) 
VIMZD^OOoOOO OO 

Kj fvj 2 ^ 1—4 l-H I H -4 M Z Z 

It || 

rgv— 4 . r— 

22ZX»hOZZ 

□ □QNNNN II uo 

^ONNNNOOO 



Q. — i UJLLJO 



O 

UJ 

J— 

< 

<y> 

a 

UJ 



O • »JO^jOZJQ^-iJO'2 OUOIS-JO 3_jO 

vt* oxx^oo-^ 11 oo<m<<m<<m< 11 <cn<r^m n<m<<m<<m 

o IS) H-h- 11 li II — ► II II II -X II -CL II -X—H ~X 1 | II — II -(Xll -CL II - 
(\||— LU — • — ^(Vj^^ir\Q^G^inQ(M^inQ^- AQ*~ *A 

II A— m— — <<A — m<f— irO- ro — - O — 

I— . LU u. GQ (XjfQ GO CD CQ OJ < UCU<t UOJ<UCU CD <UCQOUCQCa<UaJ<U OJ < 

a. I— »— h-|— t— I— »— UJ<I— U<H-UJ<i|— h-UJ<f— H-t— I— LU<h- UJ<I— OJ 

^ ^^^Mixxa:ctoo^QCUG:xu^^^o^'- < «xccoc 3 £:a:oa:a: 



in 

r- 



nT cn 



OOO OOOO 



137 



OOOOOOOOOOOOOOOOOOOOCOOOOOOOOOOOO 

flOO'O^C'O'O'O'O'ONO'O'OOOOOOOOoO^H^^^.-iHf-tHHNfM 
cm (\j cm cm <m eg c\j cm rvj c\j (\j cn c*> co on m co cn ro m co cr\ m co m rn m cn m m m m m 
OOOOOOOOOOOOOOOOQOOOOOOOOOOOOOOOO 

ooooooooooooooooooooooooooooooooo 

<<<<<<<<<<<<<« 



a 

UJ 

UJ 

^ — 
- CM 

a II 
-J 2 

a uj 

— Q 

II *» 

a. <m 

OO — * 

— ir\ 
a ii 







-J 


-J 


-J 


-J 


3 


3 


3 


3 


—4 S| 






< 


< 


< 


< 


< 


< 


< 


< 


CM 


— * 


— * 


o 


CJ 


CJ 


3 


CJ 


CJ 


CJ 


3 


3 CO 


CD 


23 


Q 


a 


a 


3 


a 


a 


3 


3 


a ^ 


3 


3 


3 


3 


3 


3 


3 


3 


3 


3 


2 3 


V 


x: 


















'D 3 


•k 


•» 


> 


O' 


o> 


'3 


0> 


O' 


O' 


JN 


II 


C3 


CO 


n 


in 


in 


in 


in 


•n 


in 


m 


a£ M 


J— 


u 


TO 


CO 


00 


00 


00 


00 


CO 


00 


U m 






—4 


—4 


—4 




>—4 


3 




H 


OO II 


•» 


«• 


I 


I 


1 


1 


1 


i 


» 


1 


II 3 


3 


-—4 


m 


in 


n 


m 


in 


n 


in 


in 


3 c 3 


> 


U. 


SJ 


M 


M 


m 


M 


M 


(M 


cM 


3 *UJ 


Sj 


■sl 


oo 


CO 


GO 


CO 


00 


X) 


oo 


oo 


2> — CC 


Si 


SI 


<— 4 


^4 


>— 4 


i-4 


—4 


-4 


-4 


-4 


-2 T3 


•k 


•k 


















sf*-* *> 


<M 


<M 


















| -CD 


UJ 


UJ 


M 


CM 


M 


<M 


<N 


M 


CM 


(M 


o -a. 


S 


2 


30 


00 


CO 


CO 


00 


00 


>00 


00 


33 !l 




*3 


















-f22< 


3 


UJ3 


CD 


CD 


CD 


CD 


CD 


3 


CD 


3 


m -cl ii 


»* 


—i - 


3 


3 


3 


31—3 


3 


3 


3 


II —4 03 


OO 


300 


< 


< 


< 


<2< 


< 


< 


<t 3 


3-wUJO 


u 


— 
















2 : 


— II X3 


Q- 


sn-o- 


00 


00 


CO 


co 2 : oo 


00 


00 


s 


2 3—00 


z 





i—4 


—4 


*3 


3 >— < 


l *— 4 


—4 






33 || > 




aao— -* 
















3 


3CUOO 




<0X00 


►1 


3 **3 


» •*! 


3 *40- 


, **3 


•■3 




. «k 


3<0 



2-U02< 

< ii <m<'0 
ex-* II -cr — 
OcM~^*nC3 i— 

< 

3 — — »<33SCL 

jaoco<^aOQ 
<ti— ^aj<tai— ^ 

CJ*— ’OCO^OULOOUJ 

o 
a 
o 
m 



cd o cdccd 
: 2 < 2 < 2 < < 



3 3 

4(X*— ‘Q- — CX n h3*-4 

• 2:>2:>2:>aj> 
< < < «-» 

<t < <ru.< 

>300300300 oo 

J —4 *JJ »— * UJ ' ' JJ —J UJ 

:azasas<s 
cj cj o 
< < <u< 



^UU^lLI^UJCL 
CD CD CDS 

:2<^<2<< 

— • — 1— — I 

-J -J JQ 

§!>X>X>d 

< < <t •— * 

< < <LL 

_J00 3o0 300 
i— II 1 j— >LU>— 1 J 

3s:3SQ2:< 

O O Cj I— 

< < CO* 



QJQO 

Q 



O 

o 

u. 

o 

(M 

3 



a 

o 

S' 



a 

s: 

3 

a 

oO 

> 

oO 

o 

3 



138 



LIST OF REFERENCES 



1. Pogue, E. W., "Experimental Observations of Geomagnetic 

Micropulsations: Land and Sea", Master's Thesis, Naval 

Postgraduate School, Monterey, September 1982. 

2. Gritzke, A. R. and Johnson, R. H., "Ocean Floor 
Geomagnetic Data Collection System", Master's Thesis, 

Naval Postgraduate School, Monterey, December 1982. 

3. Fisher, J. T., "Coherence in the Horizontal Plane of 

the Geomagnetic Field in the Frequency Range of 0-10 Hz" , 
Master's Thesis, Naval Postgraduate School, Monterey, 
December 1982. 

4. Schweiger, J. M. , "Evaluation of Geomagnetic Activity 
in the MAD Frequency Band (.04 to .06 Hz)", Master's 
Thesis, Naval Postgraduate School, Monterey, September 
1982. 

5. Jacobs, J. A., "Geomagnetic Micropulsations", p. 35-39, 
Springer-Verlay , 1970. 

6. Jackson, J. D., "Classical Electrodynamics", p. 145-148, 

John Wiley and Sons, Inc., 1962. 

7. Moose, P. A. and Thomas M. E. , "Polarization Characteristics 
of Geomagnetic Fluctuations", Transactions American 
Geophysical Union, vol . 64, no. 18, p. 212, May 1983. 



139 



INITIAL DISTRIBUTION LIST 



No. Copies 
2 



1. Defense Technical Information Center 
Cameron Station 
Alexandria, Virginia 22314 

2. Library, Code 0142 2 

Naval Postgraduate School 

Monterey, California 93940 

3. Dr. Otto Heinz, Code 61Hz 2 

Department of Physics 

Naval Postgraduate School 
Monterey, California 93940 

4. Dr. Andrew R. Ochadlick Jr., Code 610C 2 

Department of Physics 

Naval Postgraduate School 
Monterey, California 93940 

5. Dr. Paul Moose, Code 62Me 1 

Department of Electrical Engineering 

Naval Postgraduate School 
Monterey, California 93940 

6. Dr. Michael Thomas 1 

Applied Physics Lab, JHU 

Submarine Technology Division 
Johns Hopkins Road 
Laurel, MD 20707 

7. CPT Kurt B. Stevens 2 

45 Tuscarora Street 

Addison, New York 14801 

8. Office of Naval Research 1 

Attn: Mr. John G. Heacock 

Code 425 GG 

800 N. Quincy Street 

Arlington, VA 22217 

9. Naval Air Development Center 1 

Attn: Mr. Edward Yannuzzi 

Code 30 

Warminster, PA 18974 



140 



1 



10. Naval Air Systems Command 
Attn: Mr. Barry Dillon 

Code AIR-340 J 
Washington, DC 20361 



141 



